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Towards a proof of AGT conjecture by methods of 

matrix models 

A.Mironov^ A.Morozov 1 and Sh.Shakirov' 



A matrix model approach to proof of the AGT relation is briefly reviewed. It starts 



from the substitution of conformal blocks by the Dotsenko-Fateev /3-ensemble averages and 
Nekrasov functions by a double deformation of the exponentiated Seiberg-Witten prepoten- 
. tial in (3 ^ 1 and g s ^ directions. Establishing the equality of these two quantities is 

a typical matrix model problem, and it presumably can be ascertained by investigation of 
integrability properties and developing an associated Harer-Zagier technique for evaluation 
of the exact resolvent. 

Oh' 

CD ' 1 Introduction 

The AGT relation pQ experimentally discovered about 1.5 years ago, brings together a number 
of different subjects some of them listed in Fig.l. In this short review we are going to explain 
them very briefly. The details and various discussions of the AGT conjecture can be found 
£NJ , in [2]-|25j. The correspondence shown in the figure, can be described at different levels: as 
representation theory, as in the first line; as a direct equality of two quantities in the second 
line: conformal blocks of two-dimensional conformal theory and Nekrasov functions of Seiberg- 
Witten (SW) theory; and as a relation between matrix models and Seiberg-Witten theory in the 
third line. The main problem with the first two options is the lack of a reasonable conceptual 
definition of the Nekrasov functions, only operational technical definitions are available so far. 
At the third level, the Nekrasov functions are considered as providing a double-deformation 
with the help of the two e-parameters of the Seiberg-Witten theory. We find this point of view 
at the AGT conjecture the most interesting and important. It involves, at the physical level, 
the correspondence between gauge theories and integrable systems, which are well known to 
stand behind both the matrix models |26j-|29j and the SW theory [30|. l3Tj 132] . In this paper 
we rely upon the standard dictionary (gauge theory o integrable system) [31] . Also we use the 
term "matrix models" in a wider sense: the eigenvalue /3-ensembles. All the definitions from 2d 
conformal theory used below can be found in [33] and [3]. 

In a few first sections we remind the definitions of the quantities that are directly used in the 
AGT correspondence: the conformal blocks and Nekrasov functions, and formulate the AGT 
conjecture. Then we explain what definition of SW theory and of matrix models we use. After 
this, we are already able to formulate the AGT correspondence as a matrix models/SW theory 
correspondence in a more concrete form (figures 5 and 6). The rest of the paper is devoted to 
the discussion of its details. 

At the end of the paper we add three Appendices. Two of them are aimed to illustrate 
the non-widely-known chapters of matrix model theory, which are heavily exploited in the main 
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Figure 1: Interrelations between various fields touched by the AGT conjecture. 

text. The first of them is the Harer-Zagier explicit formula for the exact (all genus) one-point 
resolvent pi, which is not directly deducible from the Virasoro constraints and, thus, from the 
AMM/EO version of topological recursion; rather it is a direct corollary of a hidden integrable 
structure of the free energy. The second one is the SW representation of the /3-ensemble free 
energy, 5F = PYli §b P 1 ^ §A P l through contour integrals of the exact one-point resolvent. 
We illustrate both subjects with the simplest example of the Gaussian phase of the Hermitian 
matrix model with /? = 1. The program outlined in the present paper, requires an extension 
of this example to the Dotsenko-Fateev (DF) /3-ensemble. The third Appendix contains a brief 
description of various limits of the AGT relation, including the highly non-trivial "pure gauge" 
and "stationary" (e2 — > 0) limits. 



2 Conformal blocks 

The 4-point conformal block is a function of the double ratio of 4 points on the Riemann 
sphere, q = ^ Zl ~^ 2 w 23 ~ Z4 | which is parameterized by four external 

*" ' {z 1 -Z3)(z 2 -z 4 ) J V (a) V (1) 

and one intermediate conformal dimensions Aj (i = 1,2,3,4), A and A 2 w/ A s { ) 
the central charge c of conformal theory [33J. The conformal block is 
completely fixed by the conformal properties and, at small enough q, 

it can be expanded into the power series: Va 1 (0) Va 4 (oo) 

Ba(Ai, A 2 ; A 3 , A 4 |g) = £ q n B^ ] (1) 
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where the first coefficients of expansion are 

= 1 

B (i) = (A + Ai - A 2 )(A + A 3 - A 4 ) 
A 2A 
(2 ) = (A + Ai - A 2 )(A + Ai - A 2 + 1)(A + A 3 - A 4 )(A + A 3 - A 4 + 1) 

4A(2A + 1) 



| [(A 2 + A!)(2A + 1) + A(A - 1) - 3(A 2 - Ai) 2 ] [(A3 + A 4 )(2A + 1) + A(A - 1) - 3(A 3 - A 4 ) 2 ] 



2(2A + 1) (^2A(8A - 5) + (2A + l)c) 



These expressions become fast quite involved. However, at the particular case of large central 
charge c — > 00 the conformal block drastically simplifies becoming the hyper geometric function 

Ba(Ai,A 2 ;A 3 ,A4|?) c ^°° ^(A + Ax-A,, A + A 3 — A 4 ; 2A| g ) (2) 
3 The origin of conformal blocks 

If one writes the conformal block as a correlator in a chiral algebra with a projector to the 
intermediate state inserted 



B A (A 1 ,A 2 ;A3,A 4 |g)= (^a i( o)^a 2 (9)^^3(1)^(00)) (3) 

its generic structure can be determined using the operator product expansion of two primary 
chiral vertices 

y *w^> = E (z ^fflU (4) 

of the Virasoro algebra 

[L m , L n ] = (m - n)L m+n + r^n{n 2 - l)8 n+mfi (5) 
The conformal dimension of the vertex is given by 

L V a = A a V a (6) 
and the descendant operators are defined by 

Va,Y = L^ki ■ ■ ■ L-k-iL-k^Va (7) 
Then, the conformal block is manifestly given by the formula 



BA(Ai,A 2 ;A 3 ,A4|g)= £ q^Byy, = Y, 9 |Y| ^A 1 A 2 (^)QA(y,^)^ 3 A 4 (^) 

\Y\=\Y'\ \Y\=\Y'\ 

= £ ^ |Y| 7A 1 A 2 A(y)Q A 1 (y,y , )7A 3 A 4 A(y / ) 
\Y\=\Y>\ 



(8) 



where the 3-vertices are manifestly given by formulas 

1A U A 2 -A( Y ) = ( l -yVa(0)V Ai (1)V A2 (oo) ) = H I A + k i A 1 -A 2 + ^2k j J (9) 

i \ j<i J 

and the Shapovalov matrix is defined as a scalar product 

Qa(Y,Y') = (L_yV A (0) L_y,V A (oo)) (10) 

4 Nekrasov functions 

The Nekrasov functions for J\f = 2 SUSY SU{2) gauge theory with four fundamental matter 
hypermultiplets (the /3-function in such a theory is equal to zero) at first levels are manifestly 
given by formulas (see [M] for generic expressions) 



1 rir=i( a + /^ 



Z [2][0] 



eie 2 2a(2a + e) 

1 ]lr=i( Q -/^) 
" ~7^~ 2 " 2a(2a-e) ] 

1 Ilr=l( a + Vr)(a + Mr + £2) 



2! eie|(ei - e 2 ) 2a(2a + e 2 )(2a + e)(2a + e + e 2 ) ' 
1 llr=i (a - Mr) (« - Mr - e 2 ) 



2[ii][. 



2! eie|(ei - e 2 ) 2a(2a - e 2 )(2a - e)(2o - e - e 2 ) 

1 ]lr=l( a + /"r)( a + Mr + £l) 



Z [0][H] 



2! efe 2 (ei - e 2 ) 2a(2a + e 1 )(2a + e)(2a + e + ei) ' 

1 ]lr=l( a ~ MrKg ~ Mr ~ fl) 

2! e?e 2 (ei - e 2 ) ' 2a(2a - e 1 )(2a - e)(2a - e - ei) ' 

_J_ ]lr=l( a + / X r)(a-Mr) 

e ? e 2 ' (4a 2 -e 2 )(4a 2 -e|) ' 

„ 1 ]lr=i( Q + Mr) (a + Mr + £2) (a + Mr + 2e 2 ) 

[3][0] " 3! eie|(ei - e 2 )(ei - 2e 2 ) ' 2a(2a + e 2 )(2a + 2e 2 )(2a + e)(2a + e + e 2 )(2a + e + 2e 2 ) ' 



Here 

e = ei + e 2 

5 An origin of Nekrasov functions 

The origin of the Nekrasov functions is the LMNS integrals [35], which are integrals over a 
regularized moduli space of instantons, that is, the space of ADHM configurations. These 
integrals are manifestly given by formulas 

k J a<b [X a b ~ £ 1 J I Xab ~ £ 2 J o=l over pairs of 

Young diagrams 

(11) 
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Figure 2: The diagram describing the multi-point conformal block of the "comb" type, which is 
involved into the AGT relation. 



Figure 3: The brane diagram AGT related to the comb conformal block in Fig. 2. 



with Xab = Xa — Xb and the polynomials 

4 

P(x) = Y[(x + fi a ) , Q(x) = (x 2 - a 2 )((x + e) 2 - a 2 ) (12) 

0=1 

depend on the vacuum expectation value a of the scalar in the gauge theory and fundamental 
hypermultiplet masses [i a . These integrals can be manifestly taken and represented as sums over 
Young diagrams (integer partitions) [34j . 

The associated brane diagrams are the tropical limits of 2d Riemann surfaces. The cor- 
respondence between the conformal blocks and the brane diagrams is drawn on figures 2 and 
3. 

Here we pictured the correspondence for the comb type diagrams. The other type of diagrams 
are obtained from this one by a transformation of S-duality, see Figj4j The elements of matrix 
of the same transformation at the level of representation theory (of Virasoro algebra) are named 
the Racah coefficients (while the conformal block is nothing but the Clebsch-Gordan coefficients, 
which depend on the parameter q, since in the Virasoro case the co-product is parameterized by 
this parameter |36j). 



6 AGT relation 

The AGT conjecture literally states that 

B A (A 1 ,A 2 ;A 3 ,A4k) = J> n #i n) = Yl 9 |y| ^y'(ai,a2;a 3 ,a 4 ) = 

\Y\ = \Y'\ 

= (1 - g)- £ gl^^45 2 > - E * inf|5 X# (13) 

Y.Y> YY> 
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Figure 4: One of the diagrams obtained from the comb diagram of Fig. 2 by an S-duality trans- 
formation. 



where Z^S, are given by the formulas of section [U 

The simplest example of the conjecture comes from the first level (linear term in q): 

R (i) _ (A + Ax - A 2 )(A + A 3 - A 4 ) _ 
" 2A " 



z m + z m] + "--—■ 2a(2a + e) - — • 2a{2a _ e) + " 



which is correct provided that 



Aj = (14) 
e 

a + - = a (15) 

e e e e 

Hi = —- + a\ + Qf2, ^2 = - + ai-«2, = + °3 + «4, Ai4 = - + «3-a4, (16) 



2aia3 

z/ = 

eie2 

At higher levels one also needs 



with 



(17) 



e _ 1+ 6fa+ag =1 _ / iy 



& = \/- e i/ e 2, 5s = \/-eie2 (19) 



7 Seiberg-Witten representation of prepotential 

We call SW prepotential a function .F(ai, . . . , a n ) of n variables defined with the help of the SW 
data: a complex curve S and an meromorphic 1-differential on it, through the SW system of 
equations: 
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aj = (f> SI 

(20) 



37 
k dai 



SU(2) 

I. SW theory associated with the SU (2) gauge model. The prepotential T Q (a) of a 



In our study it appears in two contexts. 
I. SW theory associated with the S 

single variable a is defined on -£ su ( 2 ) of genus one and [371 (1-1)] 

Z su ^ : P 2 ~P+ (p)e^ - P_ (p)e _< * = -E 

where 



(21) 



+ ^ ; (1 + q) ' ^ ; + ' 2(g + l) 

and g is a function of the coupling constant r, g = 0f/#2 ( see [33 Ej). This SW system 
corresponds to the XXX spin chain on 2 sites |31j . The limit of all four hypermultiplet masses 
infinite (see s ITil) leads to the pure gauge theory on the physical side [38J, to the periodic Toda 
chain on 2 sites on the integrable side [31] and to a degenerated conformal block on the conformal 
side [4]. In this case, the corresponding SW pair (|2T1) is 



nf = pd 



(23) 



The prepotential T (a) is lifted up to the Nekrasov function J 7 ^^ (a) so that^ Sf/ ( 2 )(a) -> 
Tq U (a) in the limit ei = — e2 — > 0. 

II. SW theory associated with the matrix model. In any eigenvalue matrix model one 
can define a resolvent and its genus expansion in powers of g s (or 1/N) [39j-[42j (see s.13 for 
accurate definitions, the sign minus in front of the SW differential is due to the sign minus in 
front of the matrix model potential, (|48p ): 

_ n (mm) = p{mm){z) = / T r— L-\ = f^gf P ^ m \z) (24) 



/ ( mm ) p=0 

All pp mm ^(z) are meromorphic 1-differentials on a spectral curve x( mm ) and the free energy 
-p{mm) _ ^oo^ g 2 P -p{mm) gatisfies the gw equationSi (j^jj. We illustrate how formulas CP 
work in the simplest case of the Gaussian matrix model in Appendix II. 

In fact, along with the SW pair (^f{ mm ) ; p( mm )\ Q ne also can associate with each matrix 

model the "quasiclassical" (genus zero) SW pair \J r Q nm \ p^ m ^. 

As any partition function, the exponential of J : \ mm ) is a r-function, satisfying bilinear Hirota 
equations [33], while j^ mm ) [ s a f ar more sophisticated "quasiclassical" r-function, satisfying a 
system of highly non-linear WDVV equations [Hj 

FjFj'Fk = FkF^Fj 
where nxn matrix (Fi) JK = y^f^ ■ 
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8 DF and Selberg matrix models 



The generic Virasoro conformal block can be considered as an analytic continuation in N\ and 
N 2 of the DF type integrals [151 US] : 

$A (Ai , A 2 ; A 3 , A 4 1 g) = 

/. e oi*(0) . . e a 2 <P(q) . . e a 3 <P(l) . . ^(oo) ■ ( T ■ e H{x) . ^ V ^ . ^(y) . ^ 



^ ,0 /j j ree fields 

JVi fg AT 2 fl Ni N 2 N\ 

^ _ ^tt dXi jj / dy k T\(x t -x 3 ) 2b2 Y\(yk-yi) 2b2 T\(x l -y k ) 2b2 - 



i=l^° k=l J ° i<j k<l i,k 

J] xf\q - x t ) 2 ^\l - Xl f^ b J] yf\q - y k ) 2 ^ b (l ~ Vk?^ 

i k 



(25) 



m=l \ i / \ fc 



Integrals of this kind are often called the /3-ensemble, (3 = b 2 being the degree of the Van-der- 
Monde determinant in the integrand. 

The l.h.s. is a free chiral field correlator, evaluated with the help of the Wick theorem, and 
the r.h.s. is a double average in two Selberg models, 

1 J ° i<j 1=1 

9<y)) AT = s nI / 9(y) U(y* - y^ b2 U wf^C 1 - y^ a3b dy k 
' N2 Jo k<i ik=i 

These Selberg integrals are explicitly evaluated and are rational functions of a and b param- 
eters, decomposable into linear factors whenever / and g are the Jack polynomials sec 
[20} [T9] for details and references. 



9 A detailed scheme 

Now we are ready to draw a more detailed scheme as compared with Fig.l. With the notions 
introduced and defined above, the new picture looks like Fig.5. Below we briefly discuss the 
meaning of the three links with question marks in the low part of the table. If clarified, together 
with the already established Dotsenko-Fateev (DF) representation of conformal blocks these 
arrows would provide one possible proof of the AGT conjecture, the one based essentially on 
the matrix model technique. The structure of such a proof is outlined on another figure, Fig. 6. 

The proof of the AGT relation between the conformal block B and the Nekrasov function 
N is here actually substituted with a proof of a much simpler identity: between the genus zero 
resolvents associated with the /3-deformed spectral curves of the (Dotsenko-Fateev)matrix and 
Seiberg-Witten models, 

DFg SWp / ofi s 

This relation (j26|) can be established by a direct check, which for (3 = 1 has been performed 
in [14] . For (3^1 one still needs to know the /3-deformed version of , which is provided 
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Figure 5: The detailed diagram of the AGT conjecture. 
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Figure 6: The scheme of the proof of the AGT conjecture. 
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by the generic procedure of /3-deformation of the spectral curves, see s|TT] below. Once (|26|) is 
established, one lifts it by two canonical operations first to the equality of the full resolvents, 



PT" = P f W0 (27) 

with the help of the (/3-deformed) topological recursion and, second, to equality of the free 
energies 

-pDFp = jrSWp (28) 

with the help of the standard Seiberg-Witten procedure, which builds up the prepotential J- 
from the Seiberg-Witten differential p\. Finally, exponentiating the free energy provides the 
conformal blocks and the Nekrasov functions. In the case of conformal blocks, this equality, 

sWexpfj^vW) (29) 



is nothing but a new Dotsenko-Fateev type representation of conformal blocks with open inte- 
gration contours introduced in |18j . In the case of Nekrasov functions this statement, 



N = exp (-F^/eies) (30) 



remains to be proved, and the possibility of giving such a proof depends very much on the choice 
of a proper definition of Nekrasov functions: for various definitions see [351 145] . 

In the remaining part of this paper we briefly comment on the main ingredients of this 
suggested proof, especially on the three canonical operations: /3-deformation, topological re- 
cursion (i.e. (^-deformation) and the SW procedure. All formulas below are given for the 
simplest 4-point spherical conformal block. Generalizations to arbitrary conformal blocks, at 
least, spherical and topic, and associated quiver gauge theories are straightforward. 



10 Topological recursion 

Switching on the string coupling g s ^ 0, i.e. the deformation of the quasiclassical 

into the full J :MAMO ; has actually a functorial description, with no reference to matrix models. 
This lifting 

oo 
p=0 

or, equivalently, 

(£,f2o) — ► r— function, 

or even a "quantization" T quas i c i — > r, is now known under the name of topological recursion, 
where "topological" refers to the "genus expansion" in powers of g s . It already has a numerous 
applications to different subjects, not explicitly related to matrix models. 

The well-publicized part of the story, the AMM/EO construction |39|-[42j consists of building 
up a hierarchy of poly-differentials (multi-resolvents) on the spectral curve E. It is a part of 
a more general theory, [41] . covering also the decomposition formulas, i.e. the construction 
of integrability preserving intertwiners (of which the simplest example are VF-representations 
of [461 147]). This approach "implicitly" refers to the Virasoro constraints, i.e. the matrix 
model Ward identities [48] . and is very general. However, it is not very practical if one needs 
to construct f2, because it involves an auxiliary construction of all unneeded, and far more 
complicated, mw/ii-resolvents at intermediate steps. 
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An alternative approach [Ml [391 EH ED [52] builds 



oo 
p=0 

directly from fio, by writing and solving the difference-differential equation for fi, which involves 
recursion over N ~ l/g s (actually, over the matrix size), and follows from integrability properties 
of r = / gl- It is explicitly known only in the case of Hermitean model, but there is nothing 
preventing one from developing a similar technique of the Harer-Zagier (HZ) recursion, for other 
matrix models. 

Topological recursion: global description. Now we describe the topological recursion a 
bit more concretely. The construction implies that the whole set of multi-resolvents is intimately 
related to the U(l) current J[z) on the spectral curve E, with prescribed singularities: usually 
they are allowed at some fixed points (punctures) on E. In this approach the Virasoro constraints 
on partition function are written as 

V- (j 2 (^)) Z = J*JZ = j K{z, z') (j 2 {zfj Z = (31) 

with a certain kernel tC(z,z'), made out of the free-field Green function on E. The current is 
also "shifted": j{z) — > j(z) + AJ'(z) and partition function Z depends on the choice of: 

• the complex curve (Riemann surface) E, 

• the Green function IC(z,z'), i.e. projection operator 

• the punctures on E and associated loop operator J(z), 

• the local coordinates in the vicinity of the punctures, 

• the involution of the curve with punctures and loop operator, 

• the shift AJ(z) on E, 

• the contour C which separates two sets of punctures. 

If contour C goes around an isolated puncture, Z is actually defined by its infinitesimal 
vicinity and depends on behavior (the type of singularity) of J(z) at this particular puncture. 
Coordinate dependence is reduced to the action of a unitary operator (Bogoliubov transform, 
and exponential of bilinear function of jf) on Z. 

If contour C is moved away from the vicinity of the puncture, it can be decomposed into 
contours encircling all other punctures: this provides relations between Z's of different types, 
associated with different punctures, these are exactly the decomposition formulas. 

Topological recursion: local description. The standard recursive loop equations |39|-[42| 
are reproduced from this global construction locally in the vicinity of a fixed point. To this end, 
one has to choose the local parameter z in this vicinity and put 

J = V' + g 2 s V (32) 

where 

OO OO 

V(z) = Y,Ck(z) w V'{z) = Y J T k v k (z) (33) 

k=0 k k=0 

and Vk(z) and (k(z) are the full sets of 1-forms on E, related by the condition 

V(z)V'(z') = B{z,z') (34) 

where B(z,z') is the Bergmann kernel, i.e. (1,1) Green function B(z,z') =< d(j)(z)d(j)(z') > on 
E. The integral in (|3ip gives rise to a multiplication map x fix; — > f2s on the space of 1-forms 
fix. For hyperelliptic curves, which are double coverings of the Riemann sphere, the contour 
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C is a finite set of contours encircling the ramification points and z is the counterpart of z on 
the other sheet. Then the kernel K is actually a differential of the form which is a ratio of 
the (1,0) Green function on E (which is the primitive of the Bergmann kernel w.r.t. the second 
argument calculated from z' to z') and the Seiberg-Witten-Dijkraaf-Vafa differential 0: 

K ^ ^ = <d<t>(z)<j>(z?)>-<d(Kz) </>(?)> (35) 



See [H] and [531 S2] for details. 

Substitution of (f33|) into (pJTj) gives: 



q 4 /v 

5* E ( v * * ^ fk W n + f S (Cfc * 0) ^r^7 + 5 E K * u ')r fc Ti + 7;Tr*B )z = 

k,n>0 k,l>0 



+ \Y,( vk * + \ Tr * B z 

k,l>0 ~ J 



(36) 

where the shift of the current AJ'(z) is in charge of the *-trace of the Bergmann kernel. 
Expanding the products of 1-forms into linear combinations of £ (no v will arise due to 
projection property of the *-product), one obtains a one-dimensional set of constraints on 
log Z = X^p>o 9s P ~ 2 J~( p } . They can be also written as recurrent relations for the multiresol- 
vents 



T S ^ 



p^ m \z 1 ,...,z m ) = V(z 1 )...V(z m )T^ 
in the following form 

^ m 

p(plm+1) (z> ^ _ > Zm) = Tr *B{., •)6 p , S m , + £ B(; z { ) * p^ (., z Ifi )+ 

1 i=i 
p i 
+ E E'° (Pl|mj+1) (*' z - 7 ) *P iP2lmi/J+1) (;I/J) + -Tnp^-^+^i 



Pi=0Jd 



2 - • - 1 \ i *) ^1 J • • • > 



They are obtained simply by acting with operators V on (|3ip and putting — 5f.\ = after- 
wards. The terms with the Bergmann kernel come from the action of V on V', action on the 
V * V term gives rise to the trace of the Bergmann kernel. The notation here is as follows: the 
bullets, • mark arguments on which the *-product acts, two points are converted into a single z. If 
both bullets are arguments in the same function, we call the corresponding product *-trace, Tr*: 
for, say, H(z 1} z 2 ) = J2 m n H mnCm(z 1 )Cn(z2) the *-trace is Tr*H(; •) = £ m n H mn (( m * („)(z). 



11 /3-deformation 

/3-ensembles differ from eigenvalue matrix models [29] by substitution of the second power of 
the Van-der- Monde determinant A(M) in the measure on the space of eigenvalues by arbitrary 
power 2/3. This leads to a straightforward, but rather sophisticated deformation of integrals 
involving only traces and determinants of the integration /3-ensemble variables. The /3-character 
calculus [541 [24] which can be used to describe this deformation, automatically provides also 
the tools to handle the eigenvalue models with one external field, like Kontsevich model [55], 
GKM [27] and BGWM [28]. Character calculus in [24] is not, however, developed enough to 
describe /3-deformations of the Kontsevich and DV phases in these models as well as of their 
integrability properties. In the context of this paper it is sufficient to know the /3-deformations 



1 In simplest case of the sphere, T,h : v'h{z) = z 2 — AS corresponding to the Hermitean one-matrix model [41] 

K{z ' z ' ] = d?7^ (^) " ^)) 
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of the spectral curve and Virasoro constraints and thus of the AMM/EO topological recursion. 
Remarkably, these are the simplest things to /3-deform, and one does need to know anything 
about the /3-character calculus. In particular, formula (|25p describes the /3-ensemble with 8 = b 2 . 
This is actually a simple deformation: it affects the spectral curve S only through rescaling of 
coefficients, and it slightly changes the AMM/EO topological recursion (in particular, for (3^1 
the ^-expansion also includes odd powers of g s , f)5Q[) : contributions from open and non-oriented 
surfaces, "halves of their closed surface doubles"). 



/3-deformation of the SW data. The SW data (f20|) for the /3-ensemble is trivially deformed 



, dai 



Pi 



-P f Pi 

' Bi 



The genus zero one-point resolvent for /3-ensembles looks like (see s.13 for the details) 

-p P l ll+ W' Poll = f(z)= (W'p 
and the spectral curve is 

W'(z) 2 -ABf{z) 



W'(z)- y /W'(z) 2 -4Bf(z) J 
Po\i(z) = — dz 



2/3 



Po\i{z) 



W'(z) 
2/3 



4/3 2 



(38) 



(39) 



(40) 



For the generic SW theory one should first represent p |i a ^ P = 1 i n this form with some W'(z) 
and f(z) and then introduce /3. For example, for the SU(2) pure gauge theory (|23p (here we 
subtracted the potential part from the genus zero resolvent, see s.13) 



SU(2) 
J 0|l 



v 7 A 2 cos z — E 

75 



dz 



A cos — 
2 



E' dz 



SU(2)p 



Po\i 



~\l { A cos - 



/3E' 
(41) 



/3-deformation of the Virasoro algebra (=loop equations). The AMM/EO topological 
recursion is easily /3-deformed. For example, the Virasoro operators for Hermitean model change 
from 

d d d 2 

k a+b=n 

to 2 

dtn+2 ^ Sifc+n jrf dt a dtb y ^ A ^ >dt n ^ 6 > 

k a+b=n 

More generally, one constructs the /3-deformation of the topological recursion completely 
along the line of the previous section, with only two ingredients deformed. First of all, one has to 
deform the spectral curve (and, certainly, the related quantities like JC(z,z') and the Bergmann 
kernel) and, then, the Sugawara construction, which is behind (|3ip . The latter deformation 
is well-known and can be read off immediately using the standard conformal matrix model 
technique [56j [13] : one has to replace ([3T]) with 

V- (j\z) +(b-\) dzJj Z = (44) 



It is yet unknown how does the /3-deformation affects the integrability-inspired HZ topological 
recursion, but one can also assume that the changes will not be too drastic. 
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12 The LMNS double deformation of SW theory 

As clear from consideration of the previous sections, one can write the full free energy of the 
/3-ensemble as 

jr(mmU = to p.rec.0 [jf"^] (45) 

where the r.h.s. involves only functorial operation of topological recursion and the quasiclassical 
(genus zero) free energy. This means that the r.h.s. is fully defined by nothing but the SW pair 
(S,O ). 

Moreover, the same operation can be applied to the arbitrary SW pair (£, Qq w ) to obtain 
a (3,g s double deformation 



jrSW,fS = top rec 



(46) 



The Nekrasov functions have several different definitions, see |35t 145] for some of them. One 
of the definitions implies that they arise as a double deformation of the (exponentiated) gauge 
theory SW prepotential by two parameters e\,€2- This double deformation can presumably be 
alternatively described as switching on (3 ^ 1 and g s ^ with the usual AGT identifications 
(|19p . (3 = b 2 . In other words, one can define the Nekrasov functions as 



log 42 2) = F SU(2) >" = top-rec, [^ Ui2U ] (47) 
where at the r.h.s. stands just the usual gauge theory SW prepotential. 

13 Identification of SW and /3-ensemble differentials Qq an d 
spectral curves 

Now in order to establish the AGT conjecture one has to check first that the SW data coincides 
in the both cases: in the SW theory and in the /3-ensemble. Then, applying the topological 
recursion, one lifts the both to the full free energies which then also coincide. 

SW data in /3-ensembles. Assume one starts from the /3-ensemble 

Z ~ fUdZi exp (-- W ( Zl )] Y[( Zi - Zj fP (48) 
J i V 9s J t<j 

Then, using the invariance under the shift of integration variables <5z, ~ z ™ +1 with n > — 1, one 
immediately obtains (the double brackets here mean integration with the /3-ensemble measure 

m) 

II z n+l _ z n+l \ 

{^9. £ 1 — ~ £ W'iz^r 1 + 9 sYl ™?j = (49) 

Summing up these identities over n with the weights £ _n_2 , one rewrites these Virasoro con- 
straints in the form of a loop equation for the one-point resolvents: 

Pplit) - w'(0pi(0 + /(£) = (i - 0)</.^P - Pfin&t) (so) 

where 

W\0-W'(zi) 



(51) 
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and the subscript c means the connected correlator. Solving the loop equation (|50p in the leading 
order in the string coupling constant g s , one obtains (compare with (|40p . the sign is chosen here 
in order to guarantee the correct behaviour at large £) 



po\i(0 = ^ y bl ) 

In fact, it is usually more convenient to define p$\\{z) with the potential part subtracted, not 
changing the definition of higher order resolvents, and from now on we use the one-point resolvent 
defined exactly in such a way. The spectral curve and the SW differential are then 



S: y 



Oo = ydz 

4/3 2 



DF integrals as /3-ensembles. Choosing now the potential 

3 



W(£) = -2eiJ3a fc log(£-? fc ), { qi , q 2 , q 3 } = {0, q, 1} (54) 

fc=l 

one reproduces the /3-ensemble (125p and obtains (after rescaling y) from (153p the following SW 
data (see US (3.10)], [23 (3.38)], [HI (46)]): 

Qq = ydz 

2 _ M! M| m 2 M| + M!+m|- m 2 (1 - g) u ( 55 ) 

S: y ~ ~ + (z-1) 2 + (z-g) 2 z(z- 1) z(z-l)(z-q) 

Here 

m± = a 2 ,A = , M± = a 3) i = - (56) 

Note that the same spectral curve can be obtained via absolutely different procedure of 
studying the differential equation satisfied by the conformal block with a degenerate field inserted 
m (2-24)], [23 (88)]. 

SW theory vs. planar limit of the DF /3-ensemble. In order to compare the SW data 
(f55j) with that of the corresponding SW theory, eq. (f2Tj) , one has to make the change of variables 
in this latter 



e * ( _ fri + Mg - 2h)e + qp + - 2JQ \ „ = (57) 

which, indeed, leads to the SW pair ([55]) upon the identification 

1+? (M+-m+) 2 qpip 2 + P3P4 /.ox 

u = is = 1 

1 — q 1 — q z 1 — g 

Thus, we prove that the SW data for the DF /3-ensemble and the SU (2) gauge theory with 
four fundamental matter hypermultiplets coincide. Then they would keep to coincide after the 
topological recursion applied, i.e. the conformal block is equal to the Nekrasov functions (defined 
by the topological recursion). 
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14 The pure gauge limit 



The most familiar is the gauge SW theory in the "pure gauge" case, when the four fundamental 
masses are taken to infinity, while q — > 0, so that A 4 = qmim2msm4 remains finite. Then the 
SW pair, the spectral curve and the differential become (|23|) . In order to see this one suffices to 
make the change of variables in (|23p 



pe'* = y 



and immediately to come to 



flf = ydz 



s (S[/(2),pg . 



A 2 



E A 2 



(59) 



(60) 



y 



— + — + — 

z z z z° 



This spectral curve is exactly the curve which describes the pure gauge theory limit on the 
conformal side. 

Note that the DF /3-ensemble in the pure gauge limit looks rather different from (I25p . This 
limit is obtained by presenting the /3-ensemble integral ([25]) as an average of the product of two 
Selberg integrals and further taking the limit of these Selberg integrals. According to [24] the 
final answer is given by 



£*(A|A) 



cxp 



lim 

Ai,A2,A3,A4— >oo, q— >0 

g (A 2 -A 1 )(A 3 -A 4 )=A 4 

' ti * dttdt- k/ 



B(A 1 ,A 2 ,A 3 ,A 4 ;A,c\q) 



(61) 



(n + \t + )Z* (n-\t ) 



4=0 



where Z* (n\t\ is the partition function of the Brezin-Gross-Witten model in the character phase, 



ZAn\tk = tT^ k /k 



1 



and 



2a 

- J 

ei 



Vol fl n 



a = a 



[dU]pe 



ei + e 2 



tr U + +trW 



(62) 
(63) 



15 Conclusion 

In this paper we presented a possible strategy for proving the AGT relation between conformal 
blocks and Nekrasov functions, based on existence of the canonical double deformation of arbi- 
trary Seiberg-Witten prepotentials with the help of the topological recursion. This recursion can 
be applied both to the gauge theory SW prepotentials and to the Dotsenko-Fateev free energies. 
Since the spectral curves and SW differentials in these two cases are known to coincide, this 
actually implies the AGT relation. 

To make this argument into a constructive identification, with all steps explicit, a better 
understanding is desirable of the HZ version of the topological recursion, and thus of integrability 
properties of non-Gaussian and /3-deformed matrix models. 
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Appendix I. On HZ topological recursion 



The celebrated Hermitian model resolvent [39] (to make a contact with the standard matrix 
model notations, in this Appendix we use the definition of the resolvent without subtracting the 
potential term and rescale y by a factor of 2) 

HM , , x-y N 2lN(x 2 + N) lliV(135z 4 + 558Nz 2 + 158iV 2 ) ,„.. 
Pn \ x ) = — 7, \~ -5 H n 1 rjT\ \- ■ ■ ■ I 64 ) 

with y 2 = x 2 — AN actually satisfies the difference equation [39] 

p N+ i - 2p N + p N ^i = p' N /N (65) 
which can be also rewritten as the differential equation [501 [52] 

-Pn + V 2 Pn -x PN + 2N = (66) 

Equation ([US]) (or, equivalently, ([55]) ) is very different from the usual topological recursion. Still, 
it completely reproduces the asymptotic series (fM|) . As explained in [391 [50] such equations follow 
from integrability of the model, in this particular case from the lowest Toda-chain equation [26] 

7 HM 7 HM 02 

^N+l^N-l _ O , 7 HM (67) 

— -^hm — ~oq ogN [ ' 

where 

/ dMe- trM ^ 2+ ^ t ^M k 

Z N ~ J dMe -tr M y2 ^ 

nxn 

is the partition function of the Gaussian Hermitian matrix model. As is also explained in [50] 
(and even earlier in [39j ) it is often very convenient to consider not the matrix model at particular 
value of N, but the generating function w.r.t. N, the "grand-ensemble" partition function 

00 

Z HM (X) = £ \ N Z» M (69) 
N=0 

and its derivatives, such as the "grand-ensemble" resolvent: 

00 

p(x\X) = * N P(X\N) (70) 

N=0 

This function is in many respects more clever than its finite iV counterpart: just like pn it is 
given by an asymptotic series (see eq.(42) in [51]), but this time the series can be explicitly 
summed into an error function |39|: 

k=0 

The function p satisfies 

A - j ( — 

which is obviously equivalent to (|65p . or 

2A 

-f + ( x 2 _ AXdx) p, _ X p + ___ = (73) 

which is obviously equivalent to ([66]) . See also [39j [52], [391 EI] f° r more details. Integrability 
inspired equations of this type are often useful in various applications of matrix models (and, 
perhaps, will be useful in applications to the AGT conjecture). To illustrate their usefulness, in 
the next Appendix we use eq. ([65]) to prove the Seiberg-Witten representation of the Gaussian 
Hermitian model free energy exactly (in all genera). 



x \ X)= f A ( 1 + A ) fc (M- 1 )" = * erf (iz 1 -^) (71) 



™ x (^-^p) = f? (72) 
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Appendix II. Seiberg-Witten representation of the Gaussian- 
model partition function 

The SW representation (|20p of the matrix model free energy is not known widely enough, despite 
its extreme conceptual importance. It is only briefly mentioned even in recent matrix model 
reviews, see, for example, [571 E9 HO] • Moreover, its general prove is not actually available at 
the moment. 

In this Appendix we demonstrate how this representation works in the simplest case of the 
Gaussian Hermitian model (with f3 = 1). In this case 



, r N 2 12 '■ N ~ 1 

J Vn 

NxN NxN 



Vn J V N 

NxN 

where Vn is the volume of the unitary group. Thus 



k=l 



N 2 N ~ l 

F N = log Z N = — log 5 + Y, lo §( fc! ) ( 75 ) 

k=l 

This free energy actually gets contributions from all genera, since in genus expansion one uses 
the variables S = gN and g. Equivalently, one can keep g = 1 and consider the genus expansion 
simply as an asymptotical expansion at large N. Convenient for this purpose is the Euler- 
McLaurin summation formula [58] 

N-l N 

£/(*)= / /(^ + |i(/(iV)-/(0))+^^[/( fc - 1 )(iV)-/( fe - 1 )(0)] (76) 

k=0 q ' k=2 

i.e. 

A £ fik) = f(N) + *f {N) + *y {N) + ?±f»»(N) + ... (77) 
fc=i 

where B^ are the Bernoulli numbers [59], J2k=o Bf.x k /k\ = x/(e x — l), Bq = 1, B\ = —1/2, £?2 = 
1/6, B 3 = 0, Bi = -1/30, . . . Applying (|7l]) with 

oo „ 

/(*) = T{z + l) = (z + 1/2) log z - z + log + £ m T 1)z 2 m -i (78) 

771=1 



one finds 



^ = (iVlog(,iV)-iV) + logV^-^ + I ^3-^ + ...= 

00 5 1 ( 79 ) 
= (Nlog(gN) - N) + log - £ ^f]^ 

9=1 



Note that only the odd powers of 1/N survive in the series at the r.h.s. To obtain (|79p one makes 
use of the bilinear identity between the Bernoulli numbers, Yli+j=k (2i)\(2j)\ ^2i^2j = ~~ik^2k, 
k > 1. ° ' 

The aim of this appendix is to reproduce this expansion from the one-point resolvent (J& 



, lArt M y(z) N 2lN(z 2 + N) 11M135Z 4 + 558A^ 2 + 158iV 2 ) 
2 y b (z) y ll (z) y 17 (z) 

where y 2 (z) = z 4 — 4N. The A-period of p\ is simple: 

_ Pl (z\N)dz = N (81) 

-yflN 
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(hereafter, we use a peculiar notation § for the integral along the contour which goes around 
the points a and b). Only the first (genus-zero) term in ([80]) contributes to the j4-period, because 
the integration contour can be taken to infinity and p^ p > ~ z~ 4p_1 as z — > oo. The B-period is 
more complicated: 

dF f f°° 

— = - <p Pl {z\N)dz = -<b Pl (z\N)dz (82) 
oa Jb JvTn 

It can even seem that the contributions of higher genera diverge, since, for example, the integral 
J^Van r 2-JlatW2 a l° n g ^e ray {—y/4N, oo) diverges when k > 0. This is, however, not the 



' AN ( z 2 -AN) k / 2 

case: integrals of p\ are actually finite for k > 1, because the integration is along the contour 
surrounding the points yAN and oo, not along a segment or a ray, thus, the contour can be 
taken away from the singularities. This works exactly in the same way as in the case of the 

integral J_^]^ ( z 2_4 N ) k / 2 ; the same phenomenon also ensures the finiteness of the corrected 
Bohr-Sommerfeld periods in |12j . 

Actually, making the change of variables = one reduces the integral to the B- 
function, 

i/yHi-o^c=^§ (83) 

2 J Q T(a + b) 



in the following way: 



IN 



(z 2 -4iV)p/2 (4iV)p/2 7 s v s; s (AN)p/ 2 r(p/2) 



To see why the integral does no£ diverge, say for a > 0, one can deform the contour into a 
dumbbell, when the circle around £ = has a small radius r. Then the contribution from 
the segment diverges as §r(l — e 27rm ), while the integral along the left circle is ^ e m< ^d(f) = 
r^ e 2ma _ Clearly, these two leading divergencies cancel each other. 
With this prescription for calculation of 5-period integrals, one finds 

30 v(z) r A n(z) A 2 / N\ 

^dz = - lim <b ^dz = (NlogN - N) + 2AHog2A + O ( — ) 

'4JV 2 A^oo/^fjv 2 2 \A / 



128JV(1 - C) 5/2 ~ 12iV 



(85) 
(86) 



21iV(z 2 + iV) / 21C 7 (16-16C + 5C 2 ) ^ 1 , a7 . 



lA' 



y 11 ^) ' Jo 524288N 3 (1 - C) 11/2 12(W 3 



llJV(135z 4 + 558iVz 2 + 158AT 2 ) ^ _ Z 1 11C 11 ( 1 7280 - 34560C + 30384( 2 - 13104C 3 + 2275f 4 ) ^ _ 1 



?/ 17 (z) Jo 1073741824^5(1 - C) 17/2 252^5 

(88) 

so that the first terms of the expansion (|79j) are reproduced. Clearly, there is an identity 

Pi(z\N)dz = — — + consti + const 2 • N (89) 

B da 

i.e. the derivative of the free energy is reproduced modulo A^-constant and A^-linear terms (they 
can be removed into redefinition/rescaling of various quantities such as g, A and Z itself). 
The simplest way to check this identity in the generic form is to use equation (|65|) : 



Pl (z\N + 1) - 2 Pl (z\N) + Pl (z\N - 1) = l^ Pl ( z \N) (90) 



This equation directly implies that 

<j> Pl (z\N + l)dz - 2 j> Pl (z\N)dz + j> Pl {z\N -l)dz = ^ j> ^ Pl {z\N)dz 



(91) 
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The .B-contour integral of the full derivative does not vanish, as one could naively expect. The 
non-vanishing contribution comes from the genus zero part of the function: 

All higher genera contributions vanish. Therefore, we proves that 

/ pi(z\N + l)dz - 2 I Pl (z\N)dz + I pi(z\N - l)dz = -— (93) 
Jb Jb Jb n 

Precisely the same equation is satisfied by the r.h.s. of ([89]) : 

N-l 

F N log k\ => F N+1 - 2F N + F N -i = log N\ - \og(N - 1)! = log N (94) 

k=i 

and 

±(F N+1 -2F N + Fb- 1 )=± (95) 

Therefore, the l.h.s. and r.h.s. of (|89|) satisfy one and the same equation. This equation has 
any linear function of N as a kernel. This completes the proof. 



Appendix III. Limiting cases of AGT relation 

Conformal blocks and Nekrasov partition functions, which are identified by the AGT relation, 
depend on numerous free parameters (external and internal dimensions, e's), and one can look 
at various limits at the boundary of these moduli spaces. These limits are often highly non- 
trivial. When translated into the language of matrix models they can be also used to interrelate 
different /3-ensembles, very much in the spirit of [41J. In this Appendix we briefly review what 
is currently known about these limiting cases, see the following Table. 



limit 


AGT 


CFT 


SW theory 


/3-ensemble 


interior of moduli space 


m 


m 




DF /3-ensemble 


\,ext — > OO 

i = 1 . . . k 


® 




[601 S3 [38] 


DF /3-ensemble if k < Nc 
Unitary /3-ensembles otherwise 


A int -)■ oo 


mm 


m 


perturbative regime 


? 


c — y oo 

(9s = -> 0) 


m 






? 


e 2 ^0 


mm 




quantized SW theory 


? 


A — > degenerate value 


MM 


[331 M 







In the very left column of the Table only the quantities that becomes zero or infinite are 
written. For instance, the record A ex t — > oo in the second row implies that all the external 
dimensions and the central charge remain finite, while the record c — > oo in the fourth row 
means that all the conformal dimensions are kept finite. 

Originally the AGT conjecture relates conformally invariant theories in 2d and in 4d@ This 
means that one deals with the M = 2 SYM theory with additional adjoint, bifundamental or 
fundamental hypermultiplets, adjusted to guarantee vanishing of the /3-function (e.g. A^4 = 1 
or Np = 2Nc for the simple SU(Nc) theories). 

2 Restriction to 4d is actually inessential: dimension can be also 3, 5 or 6. 
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The conformal invariance is broken through the dimensional transmutation when masses of 
some hypermultiplets are taken to infinity. 

Another aspect of the AGT conjecture is that it occurs for the e-deformed SYM theories, 
where the Ad Lorentz invariance is violated by peculiar graviphoton backgrounds. The limit of 
9s = V — e i £ 2 ~~ > 0; leading to the ordinary SW theory, is singular. Remarkably, it is not always 
a naive planar limit in the matrix-model language: of interest are rather peculiar double-scaling 
limits, when the DF /3-ensemble reduces to another, still non-trivial /^-ensemble, not just to a 
quasiclassical approximation. 

In particular, in the pure gauge limit what occurs are unitary /3-ensembles |24j . which are 
very interesting and deserve further investigation. 

A special attention is recently attracted to the double-scaling limit €2 — > [7j. It is inter- 
esting, because in this limit a non-trivial dependence on the single e-parameter, e\, survives, 
despite g s — > 0, and the SW theory is quantized in the most naive way: by switching from 
the classical integrable system [31 J to its direct quantum counterpart, and from the spectral 
curves to the associated Baxter equation. The SW representation (|20p for the free energy now 
involves exact (quantum) Bohr-Sommerfeld periods [12]: the monodromies of the wave functions 
of (the Fourier transform of) the Baxter equation, which are simultaneously non-perturbatively 
corrected Harish-Chandra functions [25]. Moreover, the wave function itself has a wonderful 
interpretation in terms of CFT: it is equal to the conformal block with additional insertion of 
the degenerate primary [63] [TTl EH [23], [25] , associated with a surface operator insertion |8j in 
the M5/SYM language. Remarkably, such insertions into the conformal block can be also made 
for €2 7^ 0, which corresponds to switching on Whitham times, [25], still, their exact role in 
the AGT relation has to be further clarified. Also an adequate matrix model description of the 
Baxter wave function is an open problem. 
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